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SUMMARY 


A recently-introduced technique for accelerating the convergence of a 
one-step Lax-Wendroff method to a steady-state solution is discussed and its 
applicability extended to the more general class of two-step Lax-Wendroff 
methods. Several two-step methods which lead to quite efficient multiple-grid 
algorithms are discussed. 

Computational results are presented using the full two-dimensional Euler 
equations for both subcritical and shocked, supercritical flows. 

Extensions and generalizations are mentioned. 


INTRODUCTION 

Lax and Wandroff (1960) described a class of difference equations, of 
second-order accuracy in both space and time, "for approximating discontinuous 
time dependent solutions with prescribed initial data of hyperbolic systems of 
nonlinear conservation laws." Subsequently, Lax and Wendroff (1964) adapted 
these difference schemes to allow the computation of approximate weak solutions 
of multi-dimensional systems of conservation laws. 

Richtmyer (1962) devised a two-step difference scheme of Lax-Wendroff 
type. In this scheme, each step is spatially symmetric. Two- 'tep Lax- 
Wendroff schemes are inherently more efficient than one-step schemes because 
they do not require the evaluation of flux-vector Jacobian matrices. 

MacCormack (1969) introduced a two-step Lax-Wcndroff difference scheme 
in which the individual steps are not spatially symmetric. This method is a 
quite efficient member of its class. Its difference molecule contains fewer 
points than a spatially-symmetr ic scheme and it only requires information 
stored at grid points coincident with the spatial location of the final 
solutioi.. 

Lax-Wendroff algorithms have found wide acceptance in fluid mechanics 
and are commonly used both for the time-accurate computations of unsteady 
flow and for the time-asymptotic solution of steady flow problems. In the 


latter case, where accurate resolution of physical transients is not required, 
the opportunity exists to exploit this diminished restriction in order to 
accelerate convergence to the steady state. 

Ni (1981) introduced a multiple-grid technique for use in accelerating 
the convergence of a one-step Lax-Wendroff method. The present work discusses 
the extension of this technique to the wider class of two-step Lax-Wendroff 
algorithms. 

Results are presented for an application of this multiple-grid technique 
to the solution of the full two-dimensional Euler equations. This is done 
mainly for ease of computation and exposition and does not imply any restrictions 
inherent iu the technique. In particular, extensions to three dimensions and 
to the solution of the Navier-Stokes equations are intended. 


BASIC ALGORITHMS 

Consider a system of conservation laws 

It ■ 

such as the Euler equations describing inviscid flow, where : 
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Here p, u, v, p and E are respectively density, velocity components in the 
X- and y- directions, pressure and total energy per unit volume. Furthermore, 
the total energy per unit volume may be expressed as 

E = p(e + -|(u^2)) 

where the specific internal energy, e, is related to the pressure and density 
by the simple gas law 

p = (y-l)pe 

with Y denoting the ratio of specific heats. 

Schemes of Lax-Wendroff type may be arrived at intuitively by using 
Taylor's theorem to write the approximation ; 

q(t+At) = q(t) + Atq^ + q^^ (2) 
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One-Step Methods 


Since we seek solutions to Eqn. (1), time derivatives may be expressed 
as space derivatives : 
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where A and B are the Jacobian matrices : 

A = 3f/3q B = 3g/Dq 

Substitutions into Eqn. (2) results in : 

q(t+At) = q(t) ~ At(f +g ) 
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Second-order accurate spatial discretization of Eqn. (3) then yields a one-step 
Lax-Wendroff method. 

For example, we may make the following finite-volume type approximations 

j+r^^i+i,j'''^i+i,j-i^ 
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If we define the "change" in q ac cell centers, such that : 
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it then follows that 
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Consistent with the above approximations and definition, we may write the 
approximation : 
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This motivates the definitions : 
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We then make the approximations : 
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If we now define the "correction" to q at grid nodes : 
6q,. . = I q(t-fAt) - q(t) 

1 ♦ J 




(6a) 


(6b) 


we may combine Eqns. (5) and (6) to yield : 


4 



( 7 ) 


Eqns. (4) and (7) constitute the one-step Lax-Wendroff method used as a basic 
integration scheme by Ni (1981). Notice that Ni's scheme may also be thought 
of as a two-step scheme with a full time-increment predictor defined by Eqn. (4) 
and a corrector defined by Eqn. (7). However, such an interpretation of Ni's 
scheme is not totally in accord with the general practice of avoiding the 
computation of Jacobian matrices in two-step schemes. 


Two-step Metho ds 

Following Richtmyer (1962) many two-step methods which are second-order 
accurate approximations of Eqn. (3) have been developed. The two-step Lax- 
Wendroff methods of Burstein, MacCormack and Lapidus have been chosen for use 
in this study of multiple-grid acceleration both because of their ubiquity 
and their varyir.g degrees of similarity to the one-step method used by Ni. 


Burstein (1966) used a full time-increment predictor scheme which may 
be written as 
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where 


Aq 
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This method centers all quantities at the point (i, j, t + At/2) by means of 
the averaging procedure contained in Eqn. (8b). 

Notice that Eqn. (8a) is identical to Eqn. (4). Furthermore, if we 
introduce approximations of the type 
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then, for the linearized constant coefficient case, Eqn. (8b) reduces to Eqn. (7). 
Thus, under these restrictions, the methods of Burstein and Ni are equivalent 
on a step-for-step basis. 


MacCormack (1969) introduced a two-step method which may be expressed 


as 
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Here, the averaging procedure in the corrector step centers the scheme, both 
spatially and temporally, at the point (i, j, t + At/2). 

Notice that, although the methods of MacCor.nack and Ni are similar, in 
that both may be interpreted as full time-increment predictor two-step Lax- 
Wendroff schemes, they are not equivalent on a step-for-step basis. 

Lapidus (1967) used a half time-increment predictor scheme to achieve 
temporal centering without resort to the averaging necessary in the methods of 
Burstein and MacCormack. This method may be written as 
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Observe that the ch^ ge as defined by Eqn. (10a) is one-half of that 
defined by Eqn. (4) and uh. t, consequently, the correction as defined by 
Eqn. (10b) is, for the lin^irized constant coefficient case, twice that 
defined by Eqn. (7). 


tlULTIPLE-GRID ACCELERATION 

Consider the one-step Lax-Wendroff scheme, defined by Eqns. (4) and (7). 

Ni (1981) gives the following heuristic interpretation to these equations : 
the first calculates the change in q occurring in a control volume during the 
increment At , the second distributes the effects of the changes occurring in 
four nearest-nerghbor control volumes to their common central nodal point where 
they are combined to form the correction to the vector of conservation variables, 
as illustrated in Fie- 1. This interpretation leads naturally to Ni's coarse 
grid scheme, which is jllustrated in Fig. 2. 

Coarse Grid S cheme 

Wp. .’■'rill to I se successi’^ely coarser grids to propagate the finest-grid 
corrections throughout the computational domain, thus accelerating convergence 
to the steady state while maintaining the accuracy determined by the finest- 
grid discretization. Giveii a basic fine grid with the number of points in 
each direction expressible as 

2^+1 


for p a natural number, let successively coarser grids be defined by successive 
deletion of every other point in each coordinate direction. Ni's coarse grid 
solution procedure then replaces the computation of coarse grid changes by 
Eqn. (4) with a restriction of the latest finest-grid correction. This restrict- 
ed finest-grid correction is then distributed according to a coarse-grld version 
of Eqn. (7) to obtain a coarse-grid correction which is then, in turn, prolonged 
to the finest grid to become the new finest-grid correction. One time-cycle 
of Ni's multiple grid scheme is composed of an application of the one-step 
Lax-Wendroff scheme on the finest grid followed by an application of the coarse 
grid solution procedure to each successively coarser grid. 


In the basic integration scheme, a change at one grid point affects only 
its nearest neighbors while, in a k-level multiple grid scheme the same change 
affects all points up to 2^~^ mesh spacings distant. Furthermore, since the 
change is always determined by information from the finest grid and simply 
propagated by the distribution formulae for coarser grids, fine grid accuracy 
is maintained. 


Application to Two-Step Methods 

As the predictor step of two-step Lax-Wendroff schemes is analogous to 
the change as computed in Eqn. (4), a straightforward application of the coarse 
grid scheme discussed above may he used to create a simple multiple grid method. 
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As a basic flnest-grld integration scheme, we use a two-step Lax-Wendroff 
method, such as one of those discussed previously. The distribution of changes 
on the coarser grids may then be carried out using Eqn. (7), exactly as de- 
scribed for Ni's coarse grid scheme. The price to bt paid for the simplicity 
of this approach is, of course, the necessity of computing Jacobian matrices 
on the coarser grids. However, because of the vastly reduced density of coarser- 
grid points, it is advantageous to pay this price. Substantiation of this claim 
will be provided subsequently. 

More sophisticated applications to two-step methods might attempt to 
eliminate the Jacobian matrix computations by adapting the coarse-grid distribu- 
tion equation to model the corrector step of the two-step integration scheme. 

A difficulty which arises when one considers such schemes may be illustrated 
as follows. 

In a one-step method, such as Ni's scheme, the implied computation of 
f is of the form : 

'V n 

f = f + Af 

n n 
= f + A Aq 

n'\- 

= A q ; by Euler's theorem on 

homogeneous functions 

However, in a two-step method : 

% OA. 

f = Aq ; by Euler s theorem on 

homogeneous; functions 

Thus, the corrector step of such a method does not serve to distribute changes 
in the sense required for its use in a coarse grid scheme. Note that invocation 
of Euler's theorem is not essential to this arguement, but is used for purposes 
of clarity. 

'V 

One possible remedy to this problem would be simply to compute f as 

'v n^ 

f = A q 

This would mean abandoning the original objective of eliminating all of the 
Jacobian matrix computations from the coarse grid scheme. However, for certain 
two-step methods the number of matrix operations would still be less than that 
required by Ni's coarse grid scheme. 

We reserve the construction of alternative coarse grid distribution schemes 
for future study and confine ourselves here to examining the use of Ni's coarse 
grid scheme for multiple-grid acceleration of the two-step methods of Bursteln, 
MacCormack and Lapldus. 
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RESULTS 


We have computed both subcritlcal and shocked, supercritical flows in a 
straight channel with a 10% half-thick circular arc airfoil mounted on its 
lower wall. The physical problem, geometry and computational grid arc similar 
to thoseused by Rizzi and Viviand (1981), Moretti (1980), and Ni (1981). The 
computational domain is Illustrated in Fig. 3. The subcritlcal case has an 
isentropic inlet Mach number of 0.5 while that for the supercritical case Is 
0.675. 


Four computer programs were written. Implementing multiple-grid versions 
of the methods of Ni, Burstein, MacCormack and Lapidus to solve the full two- 
dimensional Euler equations as given in Eqn. (1). In each case the initial 
state, boundary conditions and (for supercritical flow) the artificial viscosity 
treatment were identical for all four programs. A sequence of four grids was 
used, as described in Table I. 


The coarse grid scheme used was as described previously and is straight- 
forward to implement for one-step methods, such as Nl's method, and for two- 
step methods using a full time-increment predictor, such as those of Burstein 
and MacCormack. Lapidus' method uses a half time-increment predictor which 
causes a slight complication, but no essential difficulty. 


To explain this complication, we consider Lapidus' method in the one- 
dimensional, linearized, constant coefficient case : 
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In the coarse grid scheme we wish to replace Eqn. (11b) with 
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In order for the coarse grid scheme to function jjroperly we tiien choose as a 
restriction operator to replace Eqn. (11a), one with the property : 


Aq = .6qr. 1 

coarse 2 fine , 

I Lapidus 


With such a choice, the coarse grid sclieme behaves as : 


Aq « 6q^. I 

coarse tine 

Ni 
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In the case of all the other algorithms, injection was used as a restric- 
tion operator. For all four algorithms, the prolongation operator used was 
linear interpolation. 

All four algorithms performed as expected, in botn the subcrltlcal and 
supercritical test cases, producing essentially Identical results. Their 
efficiencies were, however, not identical. With efficiency measured by the 
computational work required to reduce a standard error measure to a specified 
tolerance, one may conclude that the multiple-grid version of each method is 
more efficient than its single-grid counterpart. The two-step methods all 
have roughly comparable efficiencies and are substantially more efficient than 
the one-step method of N1 . This is so in spite of the fact that, on a single 
grid, the two-step methods require approximately half the computation time per 
grid point per time cycle of the one-step method. This makes application of 
Ni's coarse grid scheme relatively expensive. 

Isomach contours for the converged solutions produced by MacCormack's 
method on grid sequences of lengths one through four for the subcrltlcal and 
supercritical test cases are shown in Figs. A and 5, respectively. The 
corresponding convergence histories are shown Iti Figs. 6 and 7. We note th.it, 
for simplicity, the artificial viscosity treatment used In the supercritical 
computations with all methods was only applied on the finest grid. Hence, 
inadequate shock resolution on the coarser grids ma> have .‘flowed the convergence 
of these cases relative to their subc rit li:a 1 counterparts. 

The work associated with the artificial viscosity treatment on the finest 
giitl biases computational comparisons in favor of the multiple-grid version of 
each method for the supercritical flow test case. Accordingly, Table II 
summarizes the work required by the multiple-grid algorithms to produce a 
converged subcrltlcal solution on grid seqiunces of length one tlirough four. 

Here convergence is assumed when the averagt' ab.solute value of the total 
correction to pu over a time-cycle does not cxceid 1 x 10“^ . Notice that, 
under the conditions of this test, the ojitlmum grid sequence length is throe. 

We further note that the multiple-grid Ni algorithm requires roughly 50% more 
work than the multiple-grid MacCorm.ick algorithm. 


CONChUSION.S 

The cu.irse grid correction sciieme introduced by Ni (1981) for use with a 
one-step Lax-Wendroff method, has !)een suicessf ul 1 y a])piied to several represent- 
ative two-step Lax-Wendroff methods. 

Each resulting multiple-grid scheme Is more efficient than its single-grid 
counterp.ir t . 

The multlpl gr id oiu‘-sttq) method req ‘ ,i"S .approximately 50% more computa- 
tional work than the multiple-grid two-stiq) i..etiinds. 

Tlie possibility of constructing alternative coarse grid scht'raes h.as been 
discussed . 



£.xten8loii of the multiple-grid acceleration technique discussed here to 
three dimensions and to the Navler-Stokes equations is planned. 
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TABLE I. - GRID DESCRIPTIONS 


GRID 

1 

2 

3 

4 

NUMBER OF 
POINTS 

65 X 17 

33 X 9 

17 X 5 

9x3 


TABLE II. - WORK REQUIREMENTS FOR 

MULTIPLE-GRID ALGORITHMS 


LENGTH OF 
GRID SEQUENCE 

1 

2 

3 

4 

MacCcrmack 

1.76 

1.16 

1.00 

1.03 

Lapidus 

1.61 

1,22 

1.05 

1.13 

Burstein 

1.74 

1.24 

1.05 

1.08 

Ni 

3.10 

1.82 

1.50 

1.55 







































FIGURE 3. - Computational Domain 








